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Abstract. Collider data can play an important role in determining the parton distri¬ 
bution functions of the nucleon. I outline a formalism which makes it possible to use 
next-to-leading order calculations in such an extraction, while minimizing the amount 
of numerical computation required. 


INTRODUCTION 

The search for new physics underlying the standard model requires a precise 
understanding of known physics. At hadron colliders, this translates into the re¬ 
quirement for a precise understanding of perturbative QCD, along with knowledge 
of the non-perturbative inputs: a s and the parton distribution functions of the 
nucleon. In recent years, our detailed knowledge of perturbative QCD — in the 
form of matrix elements [1], techniques for combining virtual and real corrections in 
next-to-leading order calculations [2,3,5,4] — has grown tremendously. The same 
is true of our knowledge of a s [6], and of the parton densities at small x [7]. 

On the other hand, issues which might have been thought settled — the parton 
distributions, and in particular the gluon distribution, at moderate x — have be¬ 
come unsettled. Investigations by the CTEQ group [8], in light of the supposed 
excess of large- Et jets claimed by CDF [9], have shown that deeply inelastic scat¬ 
tering does not constrain the gluon distribution well enough to allow a claim of 
new physics to be made on the basis of the inclusive jet distribution. (The claim of 
an excess is not, however, supported by the D0 data [10].) More recently, CTEQ 
has pointed out [11] that similar statements apply to the high-x valence quark 
distributions and resulting implications of ZEUS liigh-Q 2 data [12]. 

To date, the parton densities have been taken solely as inputs in calculations of 
collider processes. To constrain the gluon distribution better, however, it would 
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be desirable to include collider data in the global fits. It is the triply-differential 
distribution [13] in clijet events, d 3 cr/dErdyid^, that will presumably give us the 
most useful information on the gluon distribution. To obtain sensible results using 
jet data, however, the use of theory at next-to-leading order is a bare minimum. 
From a theoretical point of view, this results from the presence of multiple scales: 
a jet cross section contains not only a hard-scattering scale, but also a scale char¬ 
acterizing the ‘size’ of the jet. The perturbative expansion is thus not merely an 
expansion in a s , but rather in a s In 2 y IR and a s lny IR , where yj R is a ratio of scales, 
as well. In order to ascertain that these potentially dangerous logarithms are un¬ 
der control, and that the perturbative calculation is reliable, we must calculate 
(at least) at next-to-leading order. We must therefore perform fits to collider data 
using NLO theoretical calculations. 

FITTING TO DATA 

The parton densities depend on both the momentum fraction x at which they 
are evaluated, and on a factorization scale Q 2 . Their evolution with changes in Q 2 
is governed by the Altarelli-Parisi equation [14], 

~ = Jo d y dzS ( x -y z ) p (V’Q 2 )f( z ’Q 2 )’ (!) 

whose kernel P can be computed perturbativcly. We may thus take the non- 
perturbative input to the parton densities to be their values at some fixed scale Q jj, 
and it is these values that we must fit to data. This is usually done by picking a 
parametrization of the form [17], 

fi(x, Ql ) = AiX~ Xi {l - x) ft (l + + Six). (2) 

For the traditional observables in deeply-inelastic scattering data (e.g. F 2 ), 

one has relatively simple analytic NLO formulae in terms of the parameters 
{Ai, A*, Pi, ■ ■ •}, and one can iteratively find a best fit by adjusting parameters, 
re-evaluating the theoretical prediction, and comparing with measurements. 

For jet distributions, such an iterative procedure would be extremely slow using 
existing NLO jet programs, since each iteration involves re-running the jet program, 
a matter of several hours if not a day even with present-day workstations. A better 
way of organizing the computation would thus be desirable. 

We may expect that a better organization should be possible, because most 
of the calculation — computation of the perturbative matrix elements, applying 
cuts, clustering, and the jet algorithm — doesn’t involve the parton densities in 
an essential way. I now give an outline of such a re-organization. It bears certain 
similiarities to a formalism proposed by Graudenz, Hampel, Vogt, and Berger [15] 
for jet production in deeply-inelastic scattering, but is fully general and free of 
certain limitations present in that approach. 



I shall explain the formalism in the context of a toy problem, leading-order 
glueball-glueball scattering in quarkless QCD. The same formalism, with an appro¬ 
priate sprinkling of indices and division into the different contributions that arise 
at NLO, carries over to NLO fitting of both D1S and hadron-hadron collider differ¬ 
ential cross sections. The n-jet cross section in glueball-glueball scattering, subject 
to experimental cuts, is given by 


& n. 


ff 


dxidx 2 j dLIPS^i&c + x 2 k G —■> {/cj}” =1 ) 
x f g <- G (xi, V 2 F ({ki}, x h2 ))f g ^ G (x 2 , h 2 f ({ki }, x h2 )) 


xa^{n 2 R {{ki},x lt2 ))a(gg -*• {ki}) J n ^ n ({ki}) 


(3) 


where LIPS stands for the Lorentz-invariant phase-space measure, f g ^c is the gluon 
distribution inside the glueball, and a stands for the usual leading-order partonic 
differential cross section with the running coupling a s set to 1. Note that the ki are 
implicitly dependent on x\ and x 2 as well. The renormalization and factorization 
scales Hr and hf — typically something like a jet E F — also depend on x\ i2 and 
the final-state momenta. The jet algorithm is represented by which evaluates 

to 1 if the original n-parton configuration yields n jets satisfying the experimental 
cuts, and 0 otherwise. 

It is easiest to write down the Mellin transform of the solutions to the evolution 
equation (1), using a universal evolution operator [16], 

f^o(Q 2 ) = E z (a s (Q 2 ), a 0 )f^ G (Ql ), (4) 

where z is the conjugate variable to x. An inverse Mellin transform gives us back 
the parton distribution itself, 

Vc(*, <3 2 ) = J- J c dz x-ERa,(Q\ a 0 )/ s V o (QS). (5) 

where the contour C lies to the right of all singularities in /. 

All the parameters (except ao) that we wish to fit are contained in fg<- G (Qo)- 
This function is independent of all integration variables except z, and thus can be 
pulled out of the numerical integrations in eqn. (3). The remaining z 1,2 contour 
integrals are to be performed during the fitting procedure, but this reduces to 
just a double sum of the gluon distribution function multiplied by precomputed 
numerical coefficients, and is vastly less time-consuming than the evaluation of the 
cross-section (3). That is, each step of a fitting procedure requires computing, 

I rc+ioo rc-\-io o 

-T-y/ dzj dz 2 f z l G {Ql)f z l G (Ql )^, (6) 

z ±7T J c—ioo J c—ioo 


where E Z1,Z2 are precomputed coefficients given by 



E* 1 ’* 2 = jf 1 dx x dx 2 J dLIPS( Xl k G + x 2 k' G -»• 22 

x£ 21 (a s (/4({/q}, Xi, 2 )), q;o)-S 22 ( q; s(/^f({^}, ^ 1 , 2 )), «o) (7) 

X (/^/?({^*}) •^1,2){&;}) </ n ^n({^}) • 

We see that E is a ‘cross section’ computed with x~ z E z replacing the usual parton 
distribution functions. 

To fit the remaining parameter, «o, the running coupling at Qo, we can precede 
as follows. Generate the coefficients E~ 1,Z2 for a set of cco around a “canonical” value 
(e.g. a s (Mz) = 0.118), and then fit using interpolation. While this approach would 
be vastly too time-consuming for a large number of parameters, it is acceptable for 
a lone parameter. 

As mentioned above, a similar decomposition works for differential cross sections 

in full QCD. 

The contour integrals in eqn. (6) must be performed numerically. We must first 
choose the contours, and then a set of points along each contour; the integral is 
then approximated, as usual, by a sum of the integrand evaluated at these points. 
The contour in eqn. (5) can be deformed freely into the left half-plane for z —» 00 , 
since all the singularities of the integrand are along the real axis; we can choose an 
efficient quadratic contour using the same ideas as described in ref. [18] for evolving 
parton distributions. 
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